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ABSTRACT 


A solution algorithm is presented which solves the problem of select- 
ing ports of entry for the case where the number of ports of entry is 
constrained to be small relative to the total number of ports. The 
algorithm initially considers all ports in the system considered as candi- 
dates for ports of entry and proceeds by eliminating from further consid- 
eration one port at a time until the required number is attained. An 
attempt is made to remain as close as possible to the unconstrained solu- 
tion by eliminating at each iteration the port that has the least effect 
on the objective function value. It is pointed out that the algorithm 
can yield non-optimal results in some cases but the solution is still 


better than most feasible solutions. 
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I. INTRODUCTION 


This thesis deals with a particular combinatorial problem that can be 
faced by any developing country which relies extensively on its maritime 
fleet as a principal means of transporting goods between its different 
regions. Specifically, the problem addressed is that of selecting ports 
of entry. 

Suppose a country P normally routes imported goods through a predeter- 
mined subset of its ports, henceforth referred to as ports of entry, for 
the purpose of complying with customs requirements imposed by the country. 
Country P presently has a few ports of entry which have recently become 
inadequate to handle the rapid increase of the demand for imported goods. 
Conceivably, country P can develop all of its ports to receive imported 
goods so that all goods destined for a certain region can be sent directly 
to the nearest port, thereby eliminating the added cost of transshipment. 
However, ports of entry are much more expensive to operate and maintain 
than the ordinary ports that handle only coastwise trade. Furthermore, 
the initial cost of port development might be high. As a consequence, 
the development and maintenance of all ports may be well beyond the re- 
Sources available. Hence, country P would like to increase its ports of 
entry to an acceptable number which can relieve the pressure on its 
present ports of entry, effectively serve the needs of the country and, 
at the same time, stay within a special budget appropriated for ports 
development. 

The problem faced by country P is that of determining new ports of 


entry. Since all goods imported by country P have to go through a port 
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of entry, the selection of the new ports of entry must consider the 
location of the existing ones so that the final set of ports of entry can 
effectively serve the remaining ones. A measure that can be used in 
determining the effectiveness of the system for a given choice of ports 
of entry is the total cost incurred in satisfying the total demand for 
imported goods. This total cost consists of the cost of operating and 
maintaining ports of entry and the cost of transshipping imported goods 
to their destination. The problem then reduces to the selection of ports 
of entry so that the total cost incurred in satisfying the total demand 
for imported goods is minimized. 

Section II presents a rigorous formulation of the problem as a trans- 
shipment problem with a constraint on the number of sources used. For 
convenience and simplification, the cost of developing a port into a port 
of entry is initially assumed constant and the same for any port selected. 
It is also assumed that the special budget for ports development is a 
linear multiple of the constant cost of developing a port. This facili- 
tates the determination of the number of ports of entry that can be devel- 
oped for the country. The formulation of the problem is modified slightly 
to allow for a solution using some existing network flow algorithms. 
Section III proposes an algorithm that can be used to solve the problem. 
The algorithm utilizes the Out of Kilter Algorithm of Ford and Fulkerson 
(Ref. 1) to determine the optimal flows, ignoring the ports development 
budget restriction. An example is also presented and explained. Section 
IV discusses the efficiency of the algorithm and the merits of the results 
obtained. An extension brought about by the relaxation of some of the 
assumptions is also presented and discussed. Section V summarizes the 


report and recommends some areas for further investigation. 
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ITI. PROBLEM FORMULATION 


Formally, the problem is to select some ports out of all possible 
ports through which importations can be made so that the following condi- 
tions are satisfied; 

1. The ports selected must include the existing ports of entry. In 
this thesis the number of ports selected has to be greater. than, or equal 
to, the present number of ports of entry. 

2. The demand for imported goods at each port is satisfied fully. 
This demand shall be taken as the total of all demands for imported goods 
of the region that is being served by the port. In order to justify this 
pooling of demands, it is assumed that all goods destined for any place 
in the country will pass through the port which serves that place. 

3. The quantity of goods imported through a port of entry can not 
exceed the capacity of the port. This capacity may be difficult to an 
termine due to the stochastic nature of shipment arrivals and the arrival 
of claims for goods. This report does not deal with this difficulty but 
rather assumes that a concrete value for pert capacity is available. 

-4, The total cost of satisfying the demands at all ports is minimized. 
This cost includes the total operating and maintenance cost of ports of 
entry and the total cost of transshipping the imported goods demanded at 
ordinary ports. It is assumed initially that cost is a linear multiple 
of the quantity of goods entered at or transshipped from ports of entry. 

9. The total cost of developing the ports selected can not exceed 


the special budget appropriated for ports development. 





A. DEFINITIONS AND ASSUMPTIONS 

Throughout this report, imported goods referred to are only those 
that fall under the category of general cargo. It will be assumed that 
handling charges and transshipment costs per unit of shipping is the same 
for all goods in this category. A general cargo can therefore be con- 
sidered as the single conmodity that is entered at and shipped from ports 
of entry. 

The context of any definition in this section is a graph G(N,A) com- 
posed of a non-empty set of nodes N and a set of directed arcs A that link 
the elements of N. A node is denoted by an integer. The ordered pair 
(i,j) denotes a directed arc incident from node i and incident to node j. 

Associated with each arc (i,j) are constants Cig: the cost incurred 
by a unit of flow across the arc; Lay a lower bound on the amount of flow 
that is sent across the arc; and May? the flow upper bound or arc capacity. 
Let Xs denote the actual flow sent across arc (i,j). 

The following definitions are made in connection with the formulation 
of the problem and the development of a solution algorithn; 

N = the set of all ports of the country. 

A = The set of all legs of all commercial shipping routes 

linking the ports of the country. 


Node o 


= the common single source node. This can be thought of as 
the international market which is the source of all goods 
imported by the country. 
Node n = the common single sink node. This represents the importers 
as a group who receive all goods imported. 
S = the set of all source arcs, i.e., those arcs that are 


incident from node o. 
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oj 


Lin 
Mo 
Mi 


Mp 


the set of all demand arcs, i.e., those arcs that are 
incident to node n. 

the set of all nodes in an expanded graph including the 
Single source and single sink nodes. 

the set of all arcs of an expanded graph including the 
sets of source and demand arcs. 

the special budget for ports development. 

the cost of developing part j into a port of entry. 

the increase in operating and maintenance cost per unit 
of import shipping through port j after it is made a port 
of entry. 

the cost per unit of shipping from port i to port j. 

the artificial cost for the demand at port i (assumed zero). 
the minimum quantity of goods entered at port j if it is 

a port of entry. 

a lower bound on the quantity of goods that can be shipped 
from port i to port j (assumed zero). 

a lower bound on the quantity demanded at port i. 

the maximum import capacity of port j. 

the maximum shipping capacity for the leg from port i to 
port j for 1 # 0 and j # n. 

the maximum quantity demand at port i. 

the set of all present ports of entry. 

fleet Xo 3 > 0, (0,j) « S}, the set of ports which are 
candidates for the final choice of ports of entry at 


iteration s. 





Sh 
_ 
cS 
atl 
i 


the number of elements of the set J. 


S 
Kid. = the complements of the sets K,J. 

r = the maximum number of ports of entry desired for the system. 
oe = the actual amount of goods imported through port j (or the 


amount of flow across arc (0,j)). 


As = the actual amount of goods shipped through the leg from port 
1 to port 3. 

Xen = the actual amount of the demand satisfied for port 1. 

Lag = the amount of flow allocated to the minimum cost route(s) 
found so far by the subroutine. 

ay = the least cost route from node o to node j at the end of 
iteration m. 

fe = the total cost increase per unit of shipping when flow 
across arc (0,j) is recirculated through route Rm), 

wm) = the necessary amount of flow that should be sent through 
the route i. 

oF = the minimum total cost of recirculating the flow across 


arc (0,j). 

The following assumptions are initially made in the formulation of 
the problem and the solution algorithm presented; 

1. Any-transshipment requirement in satisfying the total demand can 
be accommodated by the existing maritime fleet of the country. 

2. All goods destined for a port will be shipped via the cheapest 
Pouce. 

3. The demands at ports of entry will be satisfied by goods directly 


imported. 





4, The capacity of a port is, not less than the demand for imported 
goods of the region that it serves. 

5. All costs are linear multiples of the quantity of goods imported 
and/or transshipped. 

6. The development cost is the same for any port selected. 

7. The special budget for ports development is a linear multiple of 
the constant cost of developing any port selected. 

An extension brought about by the relaxation of some of the above 


assumptions will be discussed in Section IV. 


B. MATHEMATICAL FORMULATION 

The sets of ports N of the country together with the set of commercial 
shipping routes A that link these ports can be represented by a graph 
G(N,A). Associated with node j e N is a capacity Moy? a demand for im- 
ported goods M. and a cost Coy: Associated with each arc (i,j) « A are 
the bounds Lis = Q, A =o, and the cost per unit of shipping, Cis 

Since the required set of ports of entry are yet to be determined, 
initially consider all the ports j « N as if these are ports of entry. 
The initial graph G(N,A) would then be of the form of the graph in Figure 1 
of Appendix A. | 

The graph can be expanded to G(N*,A~) by the addition of a common 
source node 0, a common sink node n, and the sets of arcs 

See Of peace: and | — Gaiam): 1=leeween=4 } 

linking the original network with the new source and sink. The resulting 
graph G(N~,A~) will then be of the form of the graph in Figure 2. 


With each arc (0,j) e S, associate a lower bound La; =M. if jeK 


jn 
or LO = 0 if j « K, an upper bound Mo3 and a cost Cay: With each arc 
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(i,n) e€ T, associate the lower bound Li = 0, an upper bound as and a 
cost Care 0. 

If the parameters associated with the arcs of G(N’,A~) are interpreted 
as defined in Section IIA, the network obtained when flow is allowed 
through the arcs of G(N~*,A~*) can appropriately model the import system of 
the country. 

n-] 

Let Q = Me . Then Q represents the maximum total demand for the 
commodity. the ihn. as stated, can be solved using the network as a 
model by establishing a minimum cost flow pattern for Q throughout the 
network which satisfies flow conservation at nodes and which uses only 
as many arcs as the number of ports of entry required. 

The number of ports of entry required is obtained indirectly from 


condition 5 on page 6 which can be stated as 


ee 


(omipes °Y SY 
where oO and D are as defined in Section JIA and 


7 = 1 if port j is a port of entry. 
OJ 


0 otherwise. 


With assumptions 6 and 7 on page 10 , the constraint reduces to 


a or Bey eee 


(o,j)eS (mies! -= © 
indicating that the number of ports of entry can not exceed r. 
The minimum cost flow allocation of Q through the network where the 
number of source arcs used is confined to be less than or equal to r can 
then be formulated into the following mixed integer programming problem 


PI: 
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Mihm eees 


Cyak . + C. et aX. 
(ovijes 9" (Gig en omen) is 
Subject to = Xo = 
fon dvies J 
<e 
i Yunis : - 
(i k)eA- ik - (G,a)eAc J) = 0 al WP 
a © : 
(PI) (ijn)eT '” a, 
254 4 
(mies = ; 
Ls. $ a s aE ¥V(i,j)eA- 
Xo3 < Za 5M 3 : ¥V(o,j)eS 
Os Za s | La; integer 


The solution to PI, aside from determining optimal flow through the 
network, will also have r or fewer arcs (0,j)ceS with Xo; > 0. The ports 
corresponding to the nodes j to which these arcs are incident to will 
then represent the cptimal choice of ports of entry. 

Consider the problem PI without the constraint on the number of | 
source arcs (0,j)eS used. This problem reduces to the following problem 


PII: 


Minimize a 


(oejyes 05 05 * Pri ; Sper inn 
Subject me ees “ 
AOL (eye Tk ; Cer EL 29) § TEN Ase ai 
asnyer = =( 
Fe age 5 V(i,g)eA’ . 


Problem PII is nothing but a transshipment problem that can be solved 
quite easily using the Out of Kilter Algorithm of Ford and Fulkerson 
(Ref. 1). 
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The optimal solution to PII can conceivably satisfy the additional 
constraint, that is, the number of arcs (0,j)eS with Xo3 > 0 does not 
exceed r. If this is the case, then PI is also solved and the ports 
corresponding to the nodes j such that N03 > 0 are the optimal choices 
for ports of entry. If, on the other hand, the number of arcs (0,j)eS 
in the optimal solution with Xo > 0 is greater than r, then the solution 
for PII is not feasible for PI. 

Consider an optimal solution to PII with a set J = {J:X,, > omo,) jeS) 
and #(J) > r. If the flow across any #(J) - r arcs (0,j)eS are recirculated 
so that only r arcs (0,j)eS remain with XO > 0, then the resulting solu- 
tion will be feasible for PI. This property will be the basis for the 


solution algorithm that is presented in the next section. 
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TIi. SOLUTION ALGORITHM 


The solution algorithm presented in this section attempts to solve 
PI by first solving PII and deriving from it a feasible solution to PI. 
It utilizes the Out of Kilter Algorithm (Ref. 1) and an appropriate min- 
imum route algorithm (see Dreyfus (Ref. 2)) in its development. 

The solution algorithm consists of a main algorithm and a subroutine. 
The main algorithm initially solves PII to obtain the minimum cost flow 
allocation of Q through the network. This solution represents the 
optimal unconstrained solution for the import system modeled by the net- 
work. It selects the number and location of ports of entry in order to 
attain the least possible cost incurred in satisfying the demand for im- 
ported goods. The set J) = {j: Xo > 0, (0,j)eS} represents the ports 
selected. From this initial solution, the main algorithm obtains a feasi- 
ble solution to PI by reducing the number of elements of the set Jj Ome 
if this is required. This is done by recirculating the flow across 
#(J, ) - rarcs in J,, one at a time. The subroutine selects the arc 
(o,j)ed, to be eliminated at each iteration by the main algorithm by 
computing the cost of recirculating the flow across each element of the 
set J. and comparing these costs. The arc that yields the least recircu- 
lation cost is the one selected by the subroutine at that iteration. The 
algorithm terminates when the number of elements of the set Je is reduced 
Zo Y. | 
A, THE ALGORITHM 

Before starting the algorithm, construct the network G(N~*,A’) for the 


system under consideration as described in Section IIB. 
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1. Use the Out of Kilter Algorithm to determine the minimum cost 
n-] 
flow allocation of Q = => M. through the network. 


2. Eliminate all arcs (0,j) for which Koj = Q. Set L. = M. for 
gb} 1=1,...,N=le Set Sat. 
3, Determine the set 


J. ae Xai => 0; (osjyest 


for the resulting network. 
te it, 
an #(J.) <r, terminate. The solution is optimal and the ports 
corresponding to the elements of the set J. are chosen as ports of entry. 
b. #(J.) > r, use the Least Recirculation Cost Subroutine to 
determine which arc (0,j*) can be eliminated. 
5. Recirculate the flow across the arc (0,j*) by increasing the 
flow across all arcs (k,l)e “ipl by on for all m=1,2,.... Eliminate 


arc (0,j*) and go back to step 3. 


Least Recirculation Cost Subroutine 
Modify the existing network from the main algorithm by adding in 
dummy arcs (i*,3°) for each arc (j,i) with Xo Sa OetOred | lc, . eau 
Bet L.-.. = er pape =) hae ere =. 
ie ~ j x. j Ox M. : Kas and C. oe C. Delete all arcs 
(i,n) and arc (n,0o). 
1. For each arc (0,j)eS with Jed AK, obtain the total recirculation 
cost C, through the following steps: 
a. Set C. = 0, log = 0 and m=}. 
b. From the resulting network, determine the least cost route 


et from node o to node j which does not contain the arc (0,j) using an 





(m) 


appropriate minimum route algorithm. Compute Ao , the cost to recircu- 


late one unit over the route where 


_ \m) 7 ea Coan Cae, 


OF (kK TRIM 


c. If Minimum (My - Xr) SK = eee et 


(k,1)eRo OJ 9d 
(m) _ _ 
Mod Oj os 
: (m) ,,(m) 
° a ° + e e 
C, old C. ne Wo 


and select another arc (0,j). If none remain, go to step 2. 


If Minimum ier Kec (yen = Yar): 
(ky 1)eR™) k] KI oj Oj 
OJ 
iyawet Wm) = Minimum.) (M4 - Xvy) 
J (ky 1T)eR 
OJ 
“or = old *03 + Wo 
: (m) ,,(m) 
= . thas. 
C. old C. ao; Wo 
(2 Leta. = old X,, + win) for all arcs (k,1)eR)"), 


kK] k] 
(3) Eliminate arc (p,q) for which 


_ ,,(m) 
M -X = Wh, 
pq ‘pq N03 


(4) Increase m to m+] and return to step 1b. 
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2. For values of c obtained above, determine the arc (o,j*) for 
which C., = Minimum C.. This chooses the arc (0,j*) as that arc which 
Jed Nk 
can be eliminated with the least increase in flow recirculation cost. Go 


to step 5 of the main algorithm. 


B. AN EXAMPLE 

Consider a hypothetical problem where it is desired to increase the 
number of ports of entry from one to three for a system of seven ports. 
Let the graphical representation G(N,A) of the system be sy: graph in 
Figure 1 of Appendix A. Let the parameter values be those outlined in 
Table 1 of Appendix B. The problems PI and PII for the system are 
formulated in Appendix C. 

The initial network used to start the algorithm is that in Figure 2 
with K = {1}. To insure that the current port of entry will be included 
among the r=3 ports selected, the lower bound for the source arc (0,1) 
1s set equal to the demand at node 1. 

Step 1 of the algorithm requires a solution to PII for a total demand 


7 
Q = SMe = 160. The solution is shown in Figure 2 with 
i=] 


J, = {1,2,4,6,7} and #(J,) = 70% 
The solution implies that the least total cost in satisfying the demand 
Q = 160 can be achieved when the ports corresponding to the elements of 
J, are made ports of entry. However, this is more than the desired 
number. If only the elements of a proper subset of Jy are selected, an 
increase in total cost will be realized. Therefore, the subset of Jy 
must be selected in such a way that the increase in total cost is minimized. 


This is guaranteed for each step by the use of the subroutine. 


i 





The results of the subroutine for tne first iteration are summarized 
in Table 2 of Appendix B. Arc (0,6) is selected to leave at this iter- 
ation (C, = 20). 

The results of step 5 of the main algorithm is a reduced network used 


to start the next iteration. This is shown in Figure 3 with 
Jo = {1,2,4,7} and #(Jd5) = 4, 


Upon application of the subroutine for the second iteration, arc 
(0,2) is selected to leave as shown by the summary of the results in 
Table 3 of Appendix B. Recirculating the flow across the arc (0,2) results 
in the network in Figure 4. The algorithm is then terminated since 
J, = {1,4,7} and #(J3) =3=r. 

The resulting solution obtained by the algorithm can easily be shown 
to be feasible for PI by checking if the additional constraints are satis- 
fied. For the solution, Za = Zo, = Z)7 = 1 and all the rest are zero in 
order to satisfy the constraint 


constraint = Loi 
(o,j)eS J 


feasible for PI. 


;° Zo iMo3 for all arcs (0,j)eS. The 
s r is then satisfied and the solution is clearly 
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The problem, as presented, can also be solved using the partitioning 
procedure proposed by Benders (Ref. 3) or by enumerating all the possible 
combinations of r ports and selecting from these the combination that 
yields the minimum total cost in satisfying the total demand for imported 
goods. These methods are, however, not recommended when the number of 
ports in the system under consideration is fairly large while r is fairly 
small due to the amount of computational work required. In the case of 
the partitioning procedure, the computational burden can conceivably be 


great due to the size and number of integer programming problems that have 





to be solved. As for the enumeration method, the different 


-(n-r)° 
combinations that have to be investigated would just be too numerous. As 
an example, consider the case of the Philippines which has 93 national 
ports with the port of Manila as the principal import port. A current 
problem is to determine other ports which can be developed to handle 
importations in order to relieve the port of Manila of the total burden. 
Assuming that two new import ports are desired (r=3), the enumeration 
method requires the investigation of 4186 different combinations. So, for 
large networks, the algorithm proposed presents a faster method of obtain- 
ing an acceptable solution. It requires at most n-r iterations and the 
Out of Kilter Algorithm is used only once. The amount of computations 

are increased only by the number of minimum route problems that have to 

be solved at each iteration. The total amount of computational work is 
still considerably less than that required by either the partitioning 


procedure or the enumeration method. 
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The solution algorithm presented divides the set of ports at each 
jteration until a feasible subset of r ports is obtained. The division 
is done by eliminating from further consideration at each iteration, that 
port whose exclusion contributes the least increase in total cost in 
satisfying the demands for imported goods when the demand at that parti- 
cular port has to be satisfied through other ports in the set currently 
under consideration. By doing this, the algorithm will then disregard 
all the feasible sets whose elements include at least one of the ports 
eliminated. The algorithm does not necessarily yield an optimal solution 
Since there can exist a better solution involving a port previously 
eliminated. This appears to be the main drawback of the algorithm. 

Consider the example problem in Section IIIB. It can be shown that 
the solution for this problem given by the algorithm is optimal in the 
cases where r is areater than or equal to three. Suppose, instead, that 
two ports of entry are desired for the system. The solution from the 
algorithm is shown in Figure 5. The ports chosen as import ports are 
those corresponding to nodes 1 and 7 at a total cost of 2390. This, 
however, is not optimal since choosing the ports corresponding to the 
nodes | and 2 results in the least total cost. The minimum cost flow 
allocation in this case is shown in Figure 6 with a total cost of 2270. 
Nevertheless, even in this case, the solution is still a good solution 
Since it is better than all the remaining possible choices. 

In the event that an optimal solution is required, the proposed 
algorithm can be used to generate a good feasible solution, which is prob- 
ably not too far from optimal, to start other algorithms yielding optimal 
solutions such as the partitioning algorithm of Benders (Ref. 3). It is 
felt that a substantial amount of computational work can be saved if 


this is done. 20 





The proposed algorithm was developed under the assumption that de- 
velopment cost is the same for any port selected. This is more often not 
true because of the different sizes and characteristics of each port. If 
the assumption is relaxed, the number of ports of entry that can be 
created with the allocated special budget can no longer be determined 


beforehand and it becomes part of the problem. The constraint 


2 Zyast is replaced by Zz qo3 Z bei b) 
(o,j)eS (on j)eS J OJ 


The algorithm can be modified slightly to solve the resulting problem by 
obtaining at the start of each iteration the total development cost 


lid Bye 8 (o,j)ed.. 
d .2 . where Z_.. = . 
0 


(o,j)ed, ae OJ otherwise. 


and comparing this with the special budget D. If the total development 
cost is less than or equal to D, then the problem is solved and the 
elements of the set J. are selected. If, on the other hand, it is greater 
than D, the subroutine is used to select the arc (0,j*) to be eliminated. 
Upon eliminating the arc (0,j*), if the total development cost still re- 
mains larger than D, a new iteration is started. Otherwise, the algorithm 
terminates. 

An extension of the problem results when it is decided beforehand that 
the system being considered should have exactly r ports of entry. In this 


case, the special budget constraint is replaced by the constraint 


Plier >O. 
>. Z_.=yY where Z.. = OJ 


(o,5)es ou otherwise . 
After solving PII, the unconstrained solution obtained can be checked to 


determine if this new constraint is satisfied. If, 


2] 





UF = Ley ben #(J, )or and the algorithm, as presented, can 


(ond jes 
be used to obtain a solution. 


(a mee Eee then #(J))<r which implies that the number of ports 
(Ong yas 


of entry which can attain the least total cost is less than the number 
required. The algorithm can be modified to solve the problem in this 
case. In determining which port jeJ, can be added to the set J), the 
cost C. of recirculating the demand a through the arc (0,j) for all 
jed, should be calculated. Then the port j* for which 


Ces = Minimum C. 


Jed, 
can be added to the set Jy. 
S. = La = r, the unconstrained solution is also optimum for the 
(COMP IES J 


problem and the elements of the set Jy are the recommended ports. 

A final extension of the problem results when the assumption of line- 
arity of costs is relaxed. It is reasonable to expect that a port of 
entry incurs some cost even when no goods are entered through it. This 
represents the fixed minimum cost necessary to maintain the personnel and 
facilities of the port. Also, the rate of increase in cost can be ex- 
pected to decrease Prot an initially large value. However, since the 
quantity of goods entered at a port of entry is normally large, it is 
not too erroneous to assume that the rate of increase is constant. Hence, 
the operating and maintenance cost of a port can be expressed as the sum 
of a fixed cost and a linear multiple of the quantity of goods entered 
at that port. 

The addition of a fixed cost changes the problem PI to the problem 


P°-I as follows: 
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= hore eae eae 


Minimize LE a, ee ee eee Cee Oe 
(osdveS °% 9% = (fo f)eS OO) Ge cee 8 
Subject to 
Xo 9 
(Gnapies J 
= Xap > os Xs = 0 tomes ls! 
(i,k)eA’ (j,i)eAr Y 
cae = -Q 
Cane 
(P-1} p= Za att 
(O,JjJeS J 
Lis < Xi < Mi fOr alee 145 yieA- 
X93 < ZA aMo 3 fOr alleio,j)eS 
Oia La s | Za 3 integer 


where bo; is the fixed cost associated with the use of a port of entry. 
Problem P“I reduces to a fixed cost transshipment problem P“II when 

the restriction on the number of source arcs used is removed. Malek- 

Zavarei and Frisch (Ref. 4) showed that a fixed cost transshipment problem 

can be converted to a fixed cost transportation problem for which some 

good approximate solution procedures are presently available (Refs. 5 

and 6). So, it is conceivable that an algorithm similar to the one pro- 


posed in this report can be developed to solve P7I. 
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V. CONCLUSIONS 


A. SUMMARY 

The problem of selecting ports of entry from a network of ports linked 
together by commercial shipping routes was developed and an algorithm to 
solve the problem was presented. The algorithm initially considers al] 
the ports as ports of entry and gradually reduces the number-to that 
number desired for the system in such a way that the increase in cost is 
minimized. 

The algorithm presented was not intended to be the only tool on which 
to base decisions. Rather, it was meant to provide the initial ground- 
work for the development of better solution procedures and, to provide 
some initial solutions to work with. The results of the algorithm were 
Found not to be optimal in some cases but are considered good starting ~ 
points from which to arrive at better solutions. The algorithm, however, 
is useful in the sense that it substantially reduces the work necessary 
to determine a good feasible soiution especially in the case where the 
number of possible ports of entry is large enough to make the enumeration 
of all possible combinations of the ports impractical and in the cases 


where solutions not far from optimal are deemed sufficient. 


B. RECOMMENDATIONS 

The author believes that an optimal sclution to the problem can be 
attained without resorting to the enumeration method and that the solution 
process presented can be developed further to yield an optimal solution 
in all cases. An attempt should be made at determining the areas where 


the algorithm fails prior to extending the algorithm. 
24 





Finally, the extension of the problem resulting from the addition 
of fixed costs should be an interesting subject for future studies. It 
is felt that a solution process similar to the one proposed in this thesis 


can be developed for that case. 
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APPENDIX B 
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Table 2 


4 (m) 
= ROUTE ay MIN(M 057% 






Table 3 
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PI: Minimize 


Subject to 


Xo 


X18 


Xog 


X39 


X48 


X59 


X6g 


X79 


X18 


+ 


a 


4 


Xo2 


12 


Xo] 


10Xp, 


4(X15 


4(Xoq 


4(Xe5 


+ Xo3 


X13 


Xo4 


X34 


X43 


X67 


60Z 
702 
402 
802 


APPENDIX C 
EXAMPLE FORMULATION 


12X95 + 10Xp4 4 12Xo, + 20X 6 + 16Xo¢6 4 ToXp9 + 
Xoq) + 5(Xpg + Xg7) + Bg + Xeq} + O(X]7 + Xp) + 
Xa) + 8 (X56 + Xpo) + 6(X., + X43) + 3(X,5 + Xoe) + 
X7g) + O(Xyg + Xog + Xgq + Xgg + Xeg + XGg + X79) 
Xog t Xos * Xo6 * X07 = 160 
ean Wem Olle 2 ste = Gee 47 a 
NOR OSS IY NS a as = 
A Sle = 6 = 
X54 ~ 34 ~ Ao4 =u 
O55 MOS ae 
‘jeu 760 7.06 a 
S16 Mahi By Sey = Say =U 
Mp og  68le 278 
folk “wks ars" Say se 

0 < Xi, < 50 0< Xj, 6 

0 < Xoq < 20 for all i,j=1,...57 

O< X39 < 10 O< ae Seal. La integer 

0< Xag <s 20 FORCAIL J=)yc5 car 

0 < Xeq < 20 

Os X69 < 10 

0 < Xog < 30 
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PII: Minimize 10Xp) + 12Xp9 + 16X93 + 12Xqq + 2O0Xoe + 16Xo¢ + 15X59 + 
A(Xyo + Xoq) + SUK 13 + Xgq) + BUKig + Xgq) + OX, + XQ) + 
4(Xoq + Xqo) + 8UXoq + Xp) + 6(X3q + Xqg) + 3(Xe7 + X76) + 
4(Xey + X76) + OXg + Xog + Xgg + Xgg + Xgg + XGg t X7Q) 
Subject to 
i G2 0S 0A 05 06 or = We 
eee “12/0 ose “pee 7a oimeminol ~ 47 ach he 
Pomme 20 aces “ee 42 Gono ee 
gh Bag 1 i ie = 0 
Mag + Ago + Xag - Xoq - 43q ~ og amy 
en Su Syne nee 75 © “G5 - 
Toma ci G7ame le 7b) | 06 au 
fog Bp © ag Peis a ease) 9 Wee iy, SY 
Er et en, 4g 5a 468 * “76 = aL 
m5). < Xo} < 100 O< X19 < 50 ON hla << 
0 < X92 s 60 O< X59 < 20 fOr alli .d= leet 
O< X03 =< 50 O< Xn9 < 10 
Os Xoa < 60 O< Xag < 20 
Os X95 s 70 O< Keg ce 720 
0< X06 <s 40 O< Xeg < 10 
O< X97 < 80 O< Xo9 < 30 
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